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ABSTRACT 

We examine the capability of a blast-wave scenario -associated to a coronal flare or 
to the expansion of CME flanks- to reproduce a chromospheric Moreton phenomenon. 
We also simulate the Moreton event of December 06, 2006 considering both the corona 
and the chromosphere. To obtain a sufficiently strong coronal shock -able to generate 
a detectable chromospheric Moreton wave- a relatively low magnetic field intensity is 
required, in comparison with the active region values. Employing reasonable coronal 
constraints, we show that a flare ignited blast-wave or the expansion of the CME flanks 
emulated as an instantaneous or a temporal piston model, respectively, are capable to 
reproduce the observations. 
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1 INTRODUCTION 

Moreton waves, a class of large-scale chromospheric distur¬ 
bances, are detected in emission in the center and blue wing 
of the Ha spectral line, whereas they appear in absorption 
in the red wing, which is interpreted as a compression and 
subsequent relaxation of the chromosphere (Uchida 1968; 
Vrsnak et al. 2002a). They propagate forming an arc-shaped 
imprint out of the flare sites, constrained in a certain angular 
span at distances as long as 500 Mm, with radial velocities 
ranging from 500 km s -1 to 2000 km s -1 (Moreton 1960; 
Moreton & Ramsey 1960; Atliay & Moreton 1961). 

Although Moreton waves are typically observed in chro¬ 
mospheric spectral lines (Ha), there is consensus that they 
are of coronal origin since their high speeds are much larger 
than the characteristic speeds in the chromosphere. Uchida 
(1968) and Uchida et al. (1973) proposed a blast-wave sce¬ 
nario where Moreton waves are interpreted as fast-mode 
MHD shocks expanding in the corona, which produce a chro¬ 
mospheric disturbance due to the shock “sweeping” over the 
chromospheric surface. Further reinforcements of this freely 
propagating large amplitude “single wave” scenario are the 
deceleration of the wavefronts, the elongation of the pertur¬ 
bations, and the decreasing amplitude of the disturbances 
(Warmuth et al. 2001, 2004a). Commonly, the sudden oscil¬ 
lation and winking of distant filaments are associated to the 
passage of the Moreton disturbances (Gilbert et al. 2008; 


Francile et al. 2013), suggesting that Moreton waves are not 
always visible at chromospheric levels. 

At other wavelengths, similar transient wave-like fea¬ 
tures were reported, e.g., for the Helium I 10830A 
line, soft-X rays and microwaves (Gilbert & Holzer 2004; 
Warmuth et al. 2005; Vrsnak et al. 2002a; Aurass et al. 
2002; Warmuth et al. 2004a). In the extreme ultraviolet 
(EUV), global propagating disturbances were firstly ob¬ 
served by the EUV Imaging Telescope (EIT) aboard the So¬ 
lar and Heliospheric Observatory (SOHO) (Thompson et al. 
1998; Delaboudiniere et al. 1995) and, some authors argue 
that EIT waves are the coronal counterpart of Moreton 
waves as they are co-spatial (Warmuth et al. 2001, 2004b). 
Since then, they have been extensively studied giving rise 
to a controversy regarding the wave or non-wave nature 
of these “EIT” (also EUV) disturbances (see White et al. 
(2014) and references therein). Recent observations give sup¬ 
port to the coexistence of a bimodality character of these 
EIT disturbances: both the wave (fast mode) and the non¬ 
wave physical mechanisms can be at work in the same event, 
although not always detected with the current instrumen¬ 
tation (Zhukov & Auchere 2004). This would lead to the 
possibility of multiple driven mechanisms, as proposed by 
Gilbert & Holzer (2004) and Zhukov (2011). Also, while 
some authors find that the EIT and Moreton waves have 
approximately the same speed, others claim that EIT waves 
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are slower than Moreton waves by a factor of two to three 
(Chen et al. 2002, 2005; Zhang et al. 2011). 

Solar flares and coronal mass ejections (CMEs) are at¬ 
mospheric explosive phenomena capable to produce large- 
amplitude coronal disturbances and shock waves leading to 
the formation of a Moreton wave. The most straightforward 
model to give account of the shock formation is a 3D piston 
mechanism. The flare ignited blast-wave model (or simple 
3D shock model) assumes that a temporal piston mechanism 
is caused by the energy release of a flare-volume expansion 
that produces an explosion-like process driven by a pres¬ 
sure pulse (Vrsnak & Oliver 2008). On the other hand, the 
CME-driven shock model proposes that the expansion of a 
CME, together with the rise of a corresponding flux rope, 
produces a combination of a piston-shock and a bow-shock, 
which generates the large-scale shock wave (Chen & Shibata 
2000; Chen et al. 2002). Temmer et al. (2009) stated that 
the observation of Moreton waves can only be reproduced by 
applying a strong and impulsive acceleration for the source 
region expansion acting in a temporal piston mechanism sce¬ 
nario. They proposed that the expansion of the flaring region 
or the lateral expansion of the CME flanks is more likely the 
driver of the Moreton wave that the upward moving CME 
front. 

The flare-CME controversy also extends to the origin of 
type II radio burst that are usually observed close in time 
and distance to the shock source. They are interpreted to 
be the emission, at the local plasma frequency, of shock ac¬ 
celerated electrons producing Langmuir waves (Knock et al. 
2001). While there is a consensus that the emission in the 
decameter or longer wavelength range is associated with 
CMEs, metric wavelengths can be due either to flare or CME 
ignited shocks (see e.g. Vrsnak & Cliver (2008)). Recent pa¬ 
pers show evidence of type II radio burst associated with 
a flare without a CME companion (Magdalenic et al. 2012; 
Su et al. 2015). 

Thus, two different views on the origin of large-scale 
coronal shock waves arise, one favoring CMEs and the other 
preferring flares. In favor of the flare model it is argued that 
the required Moreton wave acceleration is larger or more im¬ 
pulsive than the usually observed values for CMEs. However, 
the large discrepancy between the great number of registered 
flare events and the relatively rare occurrence of observable 
coronal shocks (Cliver et al. 1999) suggests that, in addition 
to the flare explosion, another mechanism could be necessary 
to produce large-scale waves, or that a very special condition 
must be accomplished for the shock formation. In fact, fol¬ 
lowing an analysis of orders of magnitude, Vrsnak & Cliver 
(2008) showed that relatively high values of the plasma pa¬ 
rameter (plasma-to-magnetic pressure ratio) f3 « 0.1 — 0.01 
are required to ignite a coronal shock wave (at least two 
orders of magnitude larger than in an active region where 
/3 « 0.0001). This could explain the rarity of Moreton waves, 
as they should be triggered at the peripheries of active re¬ 
gions where the magnetic field is decaying. An alternative 
reasoning of this argument is given in the Appendix. 

Several numerical simulations have been carried out to 
try to explain large-scale wave formations in the solar at¬ 
mosphere. 3D MHD numerical simulations were performed 
considering a solar flare-induced pressure pulse (Wu et al. 
2001), and although the main characteristics of the ob¬ 
served EIT waves were reproduced, the plasma parameter 



Figure 1. Solid line: 2D averaged chromospheric distance ( d) 
from Q o with a 1 sigma dispersion value. Earliest wavefronts 
F\ , F 2 , F 3 before /,= 1 8:45:40 UT. A : partial quadratic fit. B : 
Total quadratic fit. C : Partial linear (it.. V) : Power-law fit (taken 
from Francile et al. (2013)). 

used /3 ~ 1, was too large. The CME scenario has also 
been simulated. In many cases an ad-hoc force is used to 
model the eruptive flux rope that triggers the expansion 
of the CME flanks and drives the shock (Chen et al. 2002, 
2005; Mei et al. 2012). The simulated shock sweeps the chro¬ 
mosphere forming the Moreton wave, where a slower wave 
(identified as the EIT wave) propagates outwards. Whether 
the resulting weak coronal CME shock can produce a de¬ 
tectable chromospheric Moreton wave is part of the contro¬ 
versy since to do so the expansion has to be accelerated to 
velocities which are rare in CMEs. 

In Francile et al. (2013) we studied a Moreton wave 
detected on December 06, 2006 with the Ha Solar Tele¬ 
scope from Argentina (HASTA) in the Ha line 656.3 nm. 
We determined the kinematics of the whole complex pro¬ 
cess through a 2D reconstruction of the HASTA and cor¬ 
responding TRACE observations (Transition Region And 
Coronal Explorer, Handy et al. (1999)). In Figure 1 we syn¬ 
thesized the observational results of the chromospheric dis¬ 
tances traveled by the Moreton wave as a 2D planar projec¬ 
tion, perpendicular to the line of sight. The evolution of the 
Moreton wave lead to the activation of two distant filaments. 
We also noted three initial irregular wavefronts that can be 
attributed to local inhomogeneities of the coronal medium 
crossed by the disturbance. We used different fits to describe 
the kinematics of the event. Curve D is a power-law fit on the 
complete set of wavefront data. The resulting initial accel¬ 
eration of the curve is ao ~ —30.2 km s~ 2 (starting with an 
initial speed of so > 2121 km s -1 at t = 18 : 42 : 30 UT). A 
partial quadratic fit trend A associated with an acceleration 
of a > —2.4 km s~ 2 was also used to adjust the data ranging 
within [18 : 45 : 40 — 18 : 49 : 02] UT with an initial speed 
of so > 1463 km s _1 (extrapolated at 18 : 42 : 28 UT). The 
quadratic acceleration value and the initial speed are higher 
than those obtained by other authors for the same event 
(Balasubramaniam et al. 2010). The partial linear fit trend 
A of the figure corresponds to a free wave speed of ~ 700 km 
s -1 . We concluded that the Moreton wave event observed on 
December 6, 2006 can be interpreted as a coronal fast-shock 
wave of a blast type originated by a single flare source dur¬ 
ing a CME ejection. However, its onset time is concurrent 
with the peak of the Lorentz force applied to the photo- 
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sphere measured by Balasubramaniam et al. (2010) showing 
an overlap with the flare explosive phase and other minor 
scale events. This argument favors the hypothesis that the 
phenomenon can be described as the chromospheric imprint 
of a fast coronal shock triggered from a single source in as¬ 
sociation with a CME. 

In this work, to reproduce the observational description 
in Francile et al. (2013), we present a 2D numerical simu¬ 
lation of the December 6 , 2006 Moreton wave. The initial 
configuration is a simplified 2D blast-wave scenario able to 
trigger a real, freely propagating MHD wave without consid¬ 
ering the magnetic field restructuring of a CME. The blast- 
wave, acting as a piston mechanism is emulated by both, 
an instantaneous pressure pulse and other one that extends 
over a short time which could also resemble the action of the 
flank expansion of a CME. The aim is to discuss whether 
the instantaneous and temporal piston models are capable 
to reproduce the main characteristics of the observations. 
We here use the word ‘instantaneous’ to distinguish a pulse 
used only as the initial condition, from another one where 
the pulse is imposed for a given lapse of time of the run, 
named ‘temporal pulse’. 


THE MODEL 

The 2D ideal MHD equations for a completely ionized hy¬ 
drogen plasma, with 7 = 5/3, (7 the ratio of specific heats) 
were implemented to study Moreton waves in the frame of 
the blast-wave scenario. We first use a pressure pulse to 
simulate the flare-volume expansion (instantaneous piston 
mechanism) that causes the blast-wave propagating fast¬ 
mode shock which sweeps the chromosphere (Uchida 1968; 
Uchida et al. 1973). The pressure pulse emulates the result 
of different complex processes that can trigger the flare ac¬ 
tivity and impulsively heat the active region (AR) (Wu et al. 
2001; Guo et al. 2013; Onofri et al. 2004). The ideal MHD 
equations, in conservative form, result: 
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where p indicates the density, v the velocity, B the magnetic 
field, p is the pressure and e the energy density given by 
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As we are interested in low coronal phenomena (flare ignition 
scenario) and considering typical large values of the coronal 
pressure scale height, « 100 Mm, we neglect the gravity 
term in the equations. 

Our aim is to describe the effect of the coronal wave over 
the transition region and the upper chromosphere. Thus, 
the chromosphere is modeled as a thin simple layer where 
the pressure is constant and the temperature and density 
abruptly change at the transition region. The height of the 
region is arbitrarily assumed as 5 Mm, approximately twice 


the height of the H Q line core formation (Leenaarts et al. 
2012). We assume that the excess of the H a core emission 
generated by the compression of the coronal shock is only 
due to the action of the upper chromosphere (observation- 
ally quantified by techniques of running differences). Fol¬ 
lowing Leenaarts et al. (2012) the emission of the H Q line 
core is strongly correlated with the mass density, and is 
only weekly modulated by the temperature and the velocity. 
Thus, the density traces the variations caused by the mag¬ 
netic field, the waves and shock waves. Also, as stated by 
Leenaarts et al. (2007) this region can be considered as op¬ 
tically thin outside dynamic magnetic structures as fibrilles. 

The shock wave characteristics are given by the coro¬ 
nal properties. Since the ratio of the coronal gas pressure to 
the magnetic pressure is /3 ~ c^/va -C 1 (c s the acoustic 
speed and va the Alfven speed) the velocity of the coronal 
fast magnetosonic shock caused by the flare will mainly de¬ 
pend on the magnetic field through the Alfven speed. Also, 
as stated by the Rankine-Hugoniot conditions, the obser¬ 
vational amplitude of the emitted wave is directly related 
to the pressure (or density) increase through the shock, i.e., 
the compression ratio (Vrsnak et al. 2002b; Vrsnak & Oliver 
2008). The question is, again, if the pressure gradient can 
produce a sufficiently strong shock wave in the ambient 
corona able to generate a detectable chromospheric More¬ 
ton wave. Although the AR pressure can increase several 
times according with the admissible values of temperature 
and density (see e.g. Aschwanden (2005)), limitations arise 
because whereas a sufficiently strong magnetic field is re¬ 
quired to reach the correct velocity of the magnetosonic 
shock, the shock intensity rapidly decays with increasing 
magnetic fields (see the Appendix). Thus, to obtain a suffi¬ 
ciently large perturbation a larger pressure pulse is needed. 
However, there is an upper limit for the temperature and 
density values. A typical AR temperature threshold should 
be as large as 40 x 10 6 K, (Aschwanden 2005). An alterna¬ 
tive would be to investigate the effect of a temporal pressure 
pulse, allowing a lower pressure pulse, though acting along 
time. With the aim of reproducing the observational HASTA 
results and considering this constraint, we analyze the influ¬ 
ence of the magnetic field in the formation of large scale 
solar waves for a uniform magnetic field assumption. 


Numerical code and initial conditions 

For the numerical simulations we use the FLASH code devel¬ 
oped at the Center for Astrophysical Thermonuclear Flashes 
(Flash Center) of the University of Chicago (Fryxell et al. 
2000). This code, currently in its fourth version, can be 
used to solve the compressible magnetohydrodynamics equa¬ 
tions with adaptive mesh refinement (AMR) capabilities. We 
choose for our simulations the “Unsplit Staggered Mesh” 
scheme (Lee et al. 2009) available in FLASH, which uses 
a high resolution finite volume method with a direction- 
ally unsplit data reconstruction and the constraint trans¬ 
port method (CT) to enforce the V ■ B = 0 condition. The 
Riemann problems of the computational interfaces are cal¬ 
culated by a Roe-type solver. 

A cartesian 2D grid with a discretization of 20 x 10 cells 
is used with 8 levels of refinement. The refinement criterion 
takes into account the variations of the density, pressure and 
magnetic field (the maximum refinement corresponds to a 
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cell of 0.39 Mm). The physical domain is set to (1000, 500) 
Mm and consists of two regions that model the solar atmo¬ 
sphere: the chromosphere formed by a small region above 
the solar surface, and the corona. 

Initially the atmosphere is in total equilibrium with 
a uniform ambient magnetic field (open-field assumption). 
Thus, the plasma pressure is constant in the initial equi¬ 
librium configuration. On the other hand, although in this 
flare-ignited scenario the interest is focused in describing 
low coronal features, the physical domain was extended in 
the vertical direction to avoid possible spurious results gen¬ 
erated by the interaction between the shock and the up¬ 
per boundary. We choose typical values of temperature and 
number density at the coronal base, T u = 1.6 x 10 6 K and 
n u = 1.2 x 10 8 cnU 3 (Wu et al. 2001). The number density 
in the chromosphere is obtained considering an initial tem¬ 
perature of Ta = 1 x 10 4 K, and the same plasma pressure as 
in the corona. As mentioned, the pressure pulse intensity is 
limited by the maximum admissible temperature and den¬ 
sity, which can increase up to T ss 40 x 10 6 K and n « 10 11 
cnU 3 e.g., if we consider a flaring loop (Aschwanden 2005). 

Firstly, the pressure pulse is instantaneously triggered 
at t = 0 s. The number density is mantained constant and 
the temperature is determined by the applied pressure pulse 
variation, but, if for a given pressure increment the maxi¬ 
mum temperature is exceeded, the density is increased to 
maintain the temperature below the threshold. The size of 
the pressure pulse is fix to 10 Mm, in accordance with typ¬ 
ical flare kernel sizes (Vrsnak & Oliver 2008). The distance 
between the pressure pulse location and the chromospheric 
surface is set to adjust the space-time interval between the 
flare event and the beginning of the Moreton wave reported 
by the observations. We assume that the flare occurs near 
the boundary of the AR where the magnetic field has already 
decay and we study the propagation of the perturbation 
from the boundary of the AR across the quiet corona where 
the magnetic field is assumed uniform. The magnetic field is 
varied within ~ [1 — 10] G to adjust the phenomenological 
values of Figure 1. Secondly, the procedure is repeated using 
a temporal pressure pulse, i.e., a pulse applied during a lapse 
of time, to be determined in order to fit the observations. 
This case could represent the action of the flank expansion 
of a CME as suggested by Temmer et al. (2009). 

Figure 2 shows the setup of the physical model, with 
the upper coronal region and the downward chromospheric 
one. We use free-flow conditions (zero gradient values) for 
the upper and lateral boundaries and a reflecting condition 
for the lower boundary to model the denser solar surface 
values. 


RESULTS AND DISCUSSION 

Kinematics: a two step mechanism 

To give account of the two step mechanism is essential to 
simulate the chromosphere as a dense thin layer where the 
coronal perturbation penetrates and rebounds. To under¬ 
stand how these two mechanisms act we analyze Figure 3. 
The figure shows a zoom, at t = 300 s, of the numerical sim¬ 
ulation domain where a coronal disturbance is triggered by 
a pressure pulse A p/p = 100 times larger than the ambient 



Figure 2. 2D simulation scheme. Corona: T u = 1.6 X 10® K, 
n u = 1.2 X 10® cm -3 , p u = 0.0265 dyn cm -2 . Chromosphere: 
Tu = 10 4 K, pu = p u , nu = 1.92x 10 10 cm -3 . The chromospheric 
height is 5 Mm and h = 35 Mm is the distance of the pressure 
pulse (Ap) from the interface corona-chromosphere. 


pressure, located at a height h = 35 Mm with a uniform 
magnetic field Bo = 1 G. There are two main effects of the 
shock front over the interface corona-chromosphere: (first 
step) an intense compression of the chromosphere acting 
persistently in the vertical direction (note the vertical dis¬ 
continuity in Figure 3a-d at x ~ 18 Mm), which is firstly ini¬ 
tiated remaining quasi-stationary, and, (second step) a cir¬ 
cular shaped shock or chromospheric disturbance, appearing 
delayed with respect to the initial compression that travels 
in the corona and “sweeps” the chromosphere. We note that 
behind the leading shock wave (x ~ 120 Mm at a coronal 
height of h ~ 35 Mm) there is a complex pattern of inter¬ 
acting waves. This nonlinear interaction is able to weaken 
the leading shock and consequently the Moreton wave (see 
a simpler study of the nonlinear interaction, where the com¬ 
bined wave front is plane in e.g., Fernandez et al. (2009), 
Costa et al. (2009) and Cecere et al. (2012)). Figure 4 shows 
a temporal scheme of the second step procedure: the initia¬ 
tion of the traveling chromospheric disturbance. At (xco, to) 
the blast coronal wave is triggered. The point (xn,ti) shows 
the space-time location where the blast arises to the inter¬ 
face corona-chromosphere (point (xci,ti) is its correspond¬ 
ing coronal counterpart). For each j, ( xcj,tj ) an d (*/j,t,) 
are a coronal and a chromospheric space-time point that 
belong to the same wavefront. Due to the geometrical char¬ 
acteristic of the scenario, while the wavefront travels from 
ti to t 2 a distance dci = xc 2 — xc l in the corona, the 
corresponding intersection, lying in the interface corona- 
chromosphere, travels a larger distance dn = X 12 — xn. 
Note that, while for earlier times dcj < dij, for longer ones 
dcj ~ dij. Therefore, the initial coronal wave speed is slower 
than the chromospheric one and they become increasingly 
similar for larger times. This corresponds to a virtual de¬ 
celeration given by the variation of the speed of a virtual 
point determined by the intersection between the chromo¬ 
sphere and the wavefront. This could give account of the 
strong initial deceleration of Moreton waves, also explain¬ 
ing why some authors find that these waves are faster than 
EIT waves (Zhang et al. 2011). Later, when the wavefront is 
capable to sweep a detectable amount of chromospheric ma¬ 
terial, a real free wave that lags with respect to its coronal 
counterpart is observed. In Figure 5a-b we display two den¬ 
sity profile slices: Figure 5a at the coronal level (h = 35 Mm 
above the interface corona-chromosphere) and Figure 5b lo¬ 
cated in the chromospheric upper layer (100 km below the 
interface corona-chromosphere), at t = 300 s. The initial 
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a) Density: p x 10 15 [cm 3 ] 


b) Pressure: p [dyn cm 2 ] 



B 



c) Velocity: v [km s x ] 

0 59 119 178 238 


i[Mm] 


d) Magnetic Field: B [G] 



0 20 40 60 80 100 120 140 160 

a: [Mm] 


Figure 3. Coronal fast MHD shock at t = 300 s with B = 1 G. 
Panel a) shows the density; panel b) shows the pressure; panel c) 
shows the velocity and panel d) shows the magnetic field. 



Figure 4. Scheme of the two step mechanism. The flare origi¬ 
nates at the corona at ( xco = 0, to)- The wavefront arrives to the 
corona-chromosphere interface at (xji = 0,ti), when the coro¬ 
nal wavefront has traveled a distance dco from the origin. At 
t 2 , the coronal wavefront has traveled a distance dci from xq i 
while the chromospheric front has traveled a distance dji > dci- 
Thus, the chromospheric wavefront speed is larger than the coro- 
nal one, dn/fa - ti) = (z /2 - xn)/(t 2 - ii) > dci /(*2 - H) = 
(xc 2 — x c\ )/ (t -2 — ii). For later times, the coronal wavefront 
speed and the chromospheric one become increasingly similar, 
dlj/{tj+1 — tj) — dcj/{tj +1 — tj). The chromospheric wavefront 
lags with respect to the coronal one. 


background magnetic field used is Bo = 3 G and the trig¬ 
gering pressure pulse is A p/p = 100. In the coronal case 
an evident shock is appreciated (200 Mm < x < 230 Mm), 
with a sharp density enhancement followed by a rarefaction. 
This fall of the density behind the shock could be an al¬ 
ternative explanation for the dimming of EUV observations 
(see also Figure 3), which are usually associated with the 
eruptive volume expansion of CMEs (Chen et al. 2002). At 



Corona 


Figure 5. Density profiles in t 
and Bq = 3 G, A p/p = 100. a) 



Chromosphere 

solar atmosphere for t = 300 s 
rona; b) chromosphere. 


the chromospheric level, we note a deep persistent vertical 
compression (first step: x < 170 Mm) that pushes down the 
interface chromosphere-corona, and a traveling front of den¬ 
sity enhancement (second step: 170 Mm < x < 220 Mm) 
more diffuse and less intense than the coronal one. Taking 
into account that the Ha opacity in the upper chromosphere 
is mainly sensitive to the mass density (the upper chromo¬ 
sphere can be considered as an optically thin media) and 
only weakly sensitive to the temperature (Leenaarts et al. 

2012 ) , we assume that this density profile gives account of 
the Moreton disturbance. As in White et al. (2014) (see Fig¬ 
ure 3 of their paper), we note that there is a lag between the 
chromospheric perturbation and the coronal signal. In this 
case, at t — 300 s, the chromospheric perturbation lags « 30 
Mm behind the coronal signal. This is reinforced by the ob¬ 
servation of the early activation of two distant filaments with 
respect to the chromospheric Moreton wave evolution, even 
in regions where it is no longer detectable, (Francile et al. 

2013) . If we consider different times (not shown in the fig¬ 
ure), after the initial transitory, we find that the two sig¬ 
nals travel with similar speeds and trajectories. In Figure 6 
we present the distance traveled by a chromospheric distur¬ 
bance as a function of time for different uniform magnetic 
field values. A typical Moreton wave kinematics is obtained 
(see Table 1), i.e., a strong initial deceleration, that gradu¬ 
ally diminishes which is larger for larger magnetosonic ve¬ 
locities (stronger magnetic fields). This characteristic behav¬ 
ior is in agreement with the results in Francile et al. (2013) 
where we detected a Moreton wave with an initial deceler¬ 
ation of ~ —30 km s -1 that diminishes to ~ —2.4 km s~ x 
in an almost quadratic trend before it finishes in a linear 
one. The t values of the table are the first detection times of 
the chromospheric traveling perturbation. The initial time 
interval with a lack of data is the time that takes the verti¬ 
cal coronal front, traveling with a fast magnetosonic shock 
speed, to arrive to the corona-chromosphere interface. How¬ 
ever, the Moreton wave detection requires that the coronal 
perturbation compresses the chromosphere beyond a certain 
threshold, which was not taken into account to perform the 
figure. As in the observations, the simulations show that the 
chromospheric wave speed gradually becomes slower until 
the shock evolves to an ordinary fast magnetosonic distur¬ 
bance (the linear region of the curves in Figure 1 and Fig¬ 
ure 6 ) (Warmuth et al. 2004b). In Francile et al. (2013) a 
compression ratio of A p > 0.08 x po was required to obtain 
a detectable chromospheric perturbation which leads to a 
certain distance between the radiant point (the probable lo¬ 
cation of the single wave source projected into the chromo- 
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Figure 6. Distances traveled by the Moreton wave as a function 
of time, for different values of the magnetic field. A p/p = 100. 


Bo[G] 

n[km s 1 ] 

a[km s 2 ] 

t[s] 

1 

594 

-7 

100 

3 

1065 

-13 

74 

5 

1876 

-35 

46 


Table 1. Simulated models: Bo is the initial background magnetic 
field, v is the velocity in x direction, a is the acceleration and t is 
the time of the first disturbance detection. 


sphere) and the place where the Moreton wave is initially 
visible. Figure 7 shows the temporal average of both, a) the 
compression ratio (Ap/po) and b) the front velocity of the 
coronal shock as a function of the pressure pulse strength 
for different magnetic field values. As mentioned before, the 
velocity and the compression ratio are mainly defined by 
the magnetic field intensity: strong fields (/3 1) cause 

large shock speeds but small compression ratios. Moreover, 
the coronal shock speed is almost constant while varying 
the pressure pulse, specially for large values of the magnetic 
field, e.g., an increase of three orders of magnitude in A p/p 
implies an enhancement of the shock speed of ~ 130% for 
a magnetic field of Bo = 1 G and, of ~ 10% for Bo = 5 
G. Note that to obtain a linear shock speed trend of about 
700 km s" 1 (as in Figure 1, which is also the typical aver¬ 
aged Moreton wave speed (Zhang et al. 2011)), we need a 
magnetic field > 3 G, but due to the rapid decrease of the 
compression ratio with increasing magnetic fields, it could be 
possible that the corresponding coronal compression ratio is 
insufficient to produce a detectable chromospheric compres¬ 
sion with the HASTA instruments (a compression ratio of 
8 % is required), even for the larger allowed values of the 
pressure pulse. This would be in accordance with the fact 
that, when detected, the Moreton events are generally asso¬ 
ciated with intense flares with an impulsive phase, usually 
leading to type II radio burst (Uchida 1968). To analyze if 
the coronal shock is capable to produce a detectable More¬ 
ton wave we study the chromospheric density profiles. We 
plot the density values at the y coordinate positions lying in 
a slice (along the x coordinate) located just below the unper¬ 
turbed upper chromospheric layer as depicted in Figure 8. 
Time t — 0 s (Figure 8a) corresponds to the flare ignition and 
subsequent times indicate the evolution of the Moreton wave 
(see the chromospheric density profiles given in Figure 8b-d, 
for t = 200, 400, 600 s, respectively). The expected morphol¬ 
ogy of the wave is reproduced, i.e., an increasing width of the 
front and a decreasing amplitude of the wave with increas- 
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Figure 7. Coronal mean compression ratio a) and velocity b) for 
the shock wave considering different values of the magnetic field 
strength and the pressure pulse intensity. 




Figure 8. Density profiles at the chromospheric surface for differ¬ 
ent times showing the evolution of the Moreton wave. The mag¬ 
netic field strength is Bo = 1 G and the pressure pulse is 50 times 
larger than the ambient pressure. 


ing time/distance. The fall in the density values behind the 
wave is related to the way the measurement is performed. As 
the density profile is obtained considering a slice just below 
the unperturbed interface corona-chromosphere, the coronal 
shock compresses the chromosphere from above enhancing 
the density (see Figure 8b) and the top of the chromosphere 
is pushed downwards. Hence, what is measured behind the 
wave is the density in the rarefaction region of the coronal 
shock, and not the density of the chromosphere. Thus, in 
order to analyze the chromospheric disturbance we evalu¬ 
ate the density profile considering an “adaptive” interface 
corona-chromosphere. This is, we adjust the measurement 
to the downward movement of the chromospheric surface. 
Also, instead of simply taking the density values consider¬ 
ing a single layer of computational cells, just below the inter¬ 
face, we construct a profile using the average density over a 
vertical region below the chromospheric surface. As the Ha 
core emission measure in the upper chromosphere is mainly 
sensitive to the density disturbance (Leenaarts et al. 2012), 
this method is better adjusted to the observations and al¬ 
lows a more precise comparison. In Figure 9a-d we show the 
density profiles for the same experiment as in Figure 8a-d, 
considering a chromospheric vertical region of d = 1 Mm 
(measured downwards from the upper layer), for the av¬ 
erage density calculation. Note, the qualitative agreement 
with the observational profiles consisting on a leading den¬ 
sity enhancement (the Moreton wave) followed by an irreg- 
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Figure 9. As Figure 8, considering an adaptive chromospheric 
surface and an average density for a vertical distance of d = 1 
Mm. The magnetic field strength is Bo = 1 G and the pressure 
pulse is 50 times larger than the ambient pressure. 


ular trend (Balasubramaniam et al. 2010) (the “activated” 
region behind the front, corresponding to the vertical shock 
compression over the chromosphere, see also the Ha line cen¬ 
ter in Figure 3 of White et al. (2014)). The differentiated two 
step behavior can be appreciated from the figure. The verti¬ 
cal line separates the static region from the other one where 
the Moreton traveling perturbation can be appreciated. 


The temporal piston 

As noted before, the coronal shock speed is mainly deter¬ 
mined by the magnetic field and almost independent of the 
pressure pulse intensity and the compression ratio (see Fig¬ 
ure 7a-b). The speeds of the chromospheric perturbations 
will also be mostly determined by the magnetic field as they 
are supposed to be the consequence of the coronal shock 
sweeping over the chromospheric surface. However, a re¬ 
quirement for the perturbations to be observed is that the 
chromosphere is compressed beyond an instrumental thresh¬ 
old. Thus, it could happen that a very strong instantaneous 
pulse, not consistent with the plasma parameters, is required 
to produce a detectable perturbation. An alternative is to 
assume that the Moreton wave is triggered by a less intense 
pulse acting for a short time, i.e., a temporal piston. 

As in Figure 7a-b, in Figure lOa-b we show the coro¬ 
nal mean compression ratio and the mean shock speed, now 
as a function of the pulse duration. Note from Figure 10b, 
that the mean shock speed is also approximately indepen¬ 
dent of the pulse duration, i.e., is fundamentally determined 
by the magnetic field value. The speed initial values of Fig¬ 
ure 10b (at t = 0 s) are the same as the ones of Figure 7b 
for an instantaneous pulse intensity of 1000. Also, varying 
the duration of the pulse, from its instantaneous value to 
another of 50 s duration the compression rate increases up 
to - 100%. 
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Figure 10. As Figure 7, coronal mean compression ratio a) and 
velocity b) for the shock wave considering different values of the 
magnetic field strength as a function of the pulse duration with 
A p/p = 1000. 

The December 6, 2006 event 

We now consider the Moreton wave event on December 6, 
2006. The aim of this analysis is to reproduce the observa¬ 
tional curve obtained in Francile et al. (2013) for the time- 
distance relation displayed in Figure 1. As seen previously, 
the main parameter required to determine the coronal shock 
speed is the magnetic field strength. To estimate this value 
we assume that the Moreton wave is due to a freely propa¬ 
gating coronal shock, which gradually decays to an ordinary 
fast-mode wave (Warmuth et al. 2001, 2004b). Considering 
the observational power law curve D of the Moreton wave 
given by Figure 1, with the wave velocity corresponding to 
the later linear trend (at t > 18:49 UT) we can calculate the 
required magnetic field strength of Bo « 3.53 G for a fast 
magnetosonic speed w/ m ~ 720 km s _1 . A distance of 35 Mm 
-between the pressure pulse location and the chromospheric 
surface- is required to reproduce the phenomenological de¬ 
lay time between the flare ignition and the emergence of 
the chromospheric disturbance, i.e., ~ 100 s. We assumed a 
coronal temperature of T = 1.6 MK and a coronal number 
density of n = 1.2 x 10 8 cm -3 . The pressure pulse inten¬ 
sity is a crucial parameter in the simulation since it needs 
to be able to compensate the effect of the large magnetic 
field (associated with lower compression ratio values, see Fig¬ 
ure 7a) and produce a sufficiently strong compression wave 
in the chromosphere to be detected by the HASTA tele¬ 
scope. The measurement of the detectable perturbation is 
performed by analyzing if the density perturbation over the 
interface corona-chromosphere is larger than the threshold 
value starting from the larger distances towards the smaller 
ones for each time. Considering these requirements we firstly 
used the maximum admisible temperature and density val¬ 
ues to set the pressure inside the flaring region which corre¬ 
sponded to A p/p ~ 10000, the larger values in Figure 7. 

Figure 11 (blue circles) shows the obtained Moreton 
wave trajectory for Bo = 3.53 G. The dashed line corre¬ 
sponds to the observational data reported in Francile et al. 
(2013) and the circles indicate the Moreton wavefront ob¬ 
tained in our simulation. We assume that radiative losses 
can be treated supposing that the upper chromospheric re¬ 
gion is an optically thin media (see e.g., Gayley & Canfield 
(1991)). Thus, we considered that the emission measure is 
proportional to the square of the particle density. The data 
indicating a Moreton front are the perturbations beyond the 
threshold. They are obtained measuring the data from the 
larger to the smaller distances. The final detection of the 
Moreton perturbation occurs when the intensity has weak¬ 
ened below the threshold (t ~ 450 s). The circles correspond- 
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Figure 11. Comparison between the observational curve and the 
numerical simulation for Bo = 3.53 G. The dashed line represents 
the 2D averaged chromospheric distance traveled by the Moreton 
wave registered on December 6, 2006 by Francile et al. (2013). 
The blue and red circles indicate the traveling wavefront obtained 
by the simulation applying an instantaneous pressure pulse of 
A p/p fs 10000 and a temporal pressure pulse of A p/p « 1000, 
applied during 40 s, respectively. A threshold of 8% in the com¬ 
pression ratio of an optically thin media is considered. 

ing to later times in Figure 11, (~ 420 s < t <~ 500 s) 
represent the stationary depression produced by the persis¬ 
tent vertical coronal compression (first step). These features 
were reported as persistent static brightenings and corre¬ 
spond to chromospheric Ha features (Delannee et al. 2007; 
White et al. 2014). In Francile et al. (2013) we find that the 
Moreton wave is detected 130 Mm away from the radiant 
point. Accordingly, in our simulations, this distance (~ 120 
Mm) corresponds to a peak of the perturbation. Considering 
the evolution at different times we obtain, as in White et al. 
(2014), that the chromospheric perturbation of December 
2006, displays a characteristic down-up vertical velocity pat¬ 
tern that lags « 50 Mm behind the coronal signal and travels 
with almost the same speed and trajectory as the coronal 
one (see Figure 5). Large values of pressure pulses can be 
expected coming from super-Alfvenic reconnection outflows 
as mentioned by Mann & Warmuth (2011). However, less 
impulsive phenomenon could be accomplished providing the 
energy to generate almost the same chromospheric perturba¬ 
tion if a less impulsive event lasts a larger time. The red cir¬ 
cles in Figure 11 correspond to the less intense temporal pis¬ 
ton case that adjusts the observational curve: A p/p « 1000, 
applied during 40 s which could resemble the action of the 
expanding flanks of a CME (Temmer et al. 2009). This pulse 
duration corresponds to the minimum value that produces 
a compression ratio beyond the instrumental threshold. 


CONCLUSIONS 

Figure 11 shows that the kinematics of the Moreton event of 
December 6, 2006 can be reproduced assuming a blast-wave 
scenario, originated from a single source, that produces a 
coronal fast mode MHD shock that sweeps the chromosphere 
generating the Moreton perturbation. This source could be 
generated by a flare that ignites the blast-wave or by the 
expansion of CME flanks due to the flare ignition. We found 
that, using typical coronal and chromospheric parameters, a 
wave of large amplitude that propagates decelerating (start¬ 
ing with an initial ao > —30 km s -2 ) with a decreasing 
amplitude until it acquires a final fast magnetosonic free 


speed linear trend of v/ m ~ 720 km s -1 is obtained. As in 
White et al. (2014) and in Francile et al. (2013) we found 
that the chromospheric perturbation displays a characteris¬ 
tic down-up vertical velocity pattern that lags behind the 
coronal signal and travels with the same speed and trajec¬ 
tory as the coronal one (see Figure 5). A uniform magnetic 
field value of Bo « 3.53 G and a distance of 35 Mm, between 
the pressure pulse location and the corona-chromospheric 
interface, are required to reproduce the phenomenological 
delay time between the flare ignition and the emergence of 
the chromospheric disturbance. The first detection of the 
Moreton wavefront is « 120 Mm far from the radiant point 
which is consistent with the value obtained in Francile et al. 
(2013). 

It has been argued that the flare explosion alone is not 
enough to give account of the Moreton event due to the dis¬ 
crepancy between the frequency of flare phenomenon and 
the relatively rare observations of coronal shocks. Thus, an¬ 
other mechanism would be necessary or a very special con¬ 
dition must be accomplished to give account of the shock 
formation that produces the large-scale Moreton event. If 
the modeling used here is an accurate one, our results show 
that, given typical coronal and chromospheric parameters, a 
set of limiting conditions are required to obtain both, a final 
linear trend associated with a fast magnetosonic speed and a 
strong pressure pulse able to generate a detectable compres¬ 
sion ratio. The final fast magnetosonic free speed requires a 
definite magnetic field intensity (Figure 7b), but if the in¬ 
tensity of the magnetic field is large the compression ratio 
values are small and hinder the wave detection (Figure 7a). 
Accordingly, as seen in the Appendix, the flare expansion 
should be located at the periphery of the AR in order to 
avoid large magnetic field values (where /3 ~ 1 x 10 -4 ). Also, 
the fact that a very strong pressure pulse (limited by the ad¬ 
missible values of temperature and density) was required to 
obtain an emission enhancement of only 8% for later times 
reinforces the argument. 
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APPENDIX A: IMPORTANCE OF THE 0 
PARAMETER TO IGNITE A CORONAL 
SHOCK 

Consider a steady shock plane-parallel to the magnetic 
field (Bo) and normal to the ^-coordinate. Making use of 
the Rankine-Hugoniot conditions, the conservation laws for 
mass, momentum, energy and magnetic flux read (see for 
example Torrilhon (2003) and Draine (2011)): 


povo = pivi, 
vqBo = viBi, 

Po v o + Po + = 2 

no 


2 | . B? 

ST - P lV 1 +P 1 + ST’ 


\povl + ^Vopo + vof^-] = 

= Vl [ipiuf + ^JVIPI + Vl 


(Al) 


where the subscripts 0 and 1 denote preshock and postshock 
conditions, respectively, v is the velocity in ^-direction in a 
frame fixed to the shock front. 

The non-trivial solution is one that produces different 
values on each side of the discontinuity. Considering the 
preshock plasma at rest, it results that vo = v s in the moving 
frame, where v s is the shock speed. Defining the compression 
ratio X = pi/po (Draine 2011): 

2(2 - 7 )X 2 + 7 [(7 - tyPMs + 2(1 + 0)] X- 

7(7 + l)/SATs = 0 


(A2) 


M a = Vs/cs is the shock Mach number. 

As the shock speed is equal or larger than the fast mag- 
netosonic speed Vf m = > a shock compression 

implies v s = avfm, where a > 1, thus: 


M 2 = a 2 ^±^ 


Tp + \ 

We can rewrite Eq. A2 in terms of the parameter 0 

2(2 - 7 )X 2 + [a 2 (7 - 1)(2 + 7 0) + 2(1 + 0)\ X- 
a 2 ( 7 + 1)(2 + 7/S) =0. 


(A3) 


(A4) 


Note that, in the case of interest (1 < 7 < 2), Eq. A4 
has only one positive root which satisfies X > 1. The limit 
X —¥ 1, is obtained when 0 —> 0, independently of the values 
of 7 and a. Thus, the flare-associated pressure pulse cannot 
ignite a shock wave in strong field regions where 0 —> 0. A 
detailed study on MHD shock wave formation can be found 
in Vrsnak & Lulic (2000a) and Vrsnak & Lulic (2000b). 

This paper has been typeset from a TpX/lXTpX file prepared by 
the author. 
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